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SUMMART 


A basic theory of generalized linearized stipersonic conical flow 
for both inside and outside the Mach cone was developed and applied to 
several specific problems including unsteady-flow conditions. A tri- 
angular lifting wing in pitching and rolling with both subsonic and 
si; 5 )ersonic leading edges was investigated and pressure coefficients 
were obtained. A family of thin sweptback triangiilar wings having sym- 
metrical thickness distribution was also investigated and analytic 
expressions for wave drag and pressure coefficients were determined. 
Values of wave drag coefficients were calculated and the results pre- 
sented graphically. This theory stems from a fundamental idea of 
Mr. G. N. Ward. 


INTRODUCTION 


Of the methods and theories of linearized supersonic flow one of 
the most productive of results directly applicable to aerodynamic prob- 
lems has been the theory of conical flow. Conical flow refers to flow 
in >diich the pressure perturbation and velocity components are constant 
along straight lines or rays passing through a fixed pointj for such 
flow the assumption of linearity is not needed and nonlinear conical 
flows are of great practical interest. In the generalization of conical 
flow presented in this report, however, the linearizing assumptions are 
necessary and the usual wave equation of steady linearized supersonic 
flow will be considered. 

In generalizing the concept of conical flow attention is focused 
on the homogeneity of the solutions' with the vertex as origin. In 
regular conical flow the solutions for the velocity conponents are 
homogeneous of degree 0, while the corresponding velocity potential 
whose gradient is the vector velocity is homogeneous of degree 1. In 
generalized conical flow solutions may be considered for which the 
velocity potential is homogeneous of degree n with the velocity com- 
ponents homogeneous of degree n - 1. The quantity n, called the 
degree of the generalized conical flow, may usually be an integer but 
does not need to be. A further generalization comes from supeiposition 
of solutions of different homogeneity. 



2 


MCA TN 2667 


It is the ptii^ose of this report to present the basic theory of 
generalized linearized conical flow as developed from a fundamental idea 
of Ward and to illustrate the theory by some aerodynamic applications. 

The particular elements of linearized supersonic-flow theory which will 
be needed are presented in brief. The point of view used in the analysis 
is that used in reference 1 , which thus may be considered a basic refer- 
ence for this report. Some of the essential steps in the development 
are the saiEe as those in the conical-flow theory of Goldstein and Ward 
(reference 2) . For a careful and thorough investigation and discussion 
of the conical-flow theory the reader is also referred to the work of 
Lagerstrom (reference 3 ) . For an investigation of generalized conical- 
flow theory from quite a different point of view reference may be made 
to the theory of Germain (reference i|.) . No attempt is made to give more 
than the most recent references as quite complete bibliographies may be 
found in the references which are cited. 

The coordinate system used is the right-handed system used in 
references 1 to 3 "Hith' the third coordinate in the free-stream direction. 
There are several advantages rtiich may be claimed for this particular 
orientation: The timelike position in the wave equation is delegated 

to the last coordinate; in accordance with customary notation with this 
equation; the variables x and y appear on an equal basis; and liith 
X taken as the lateral variable a flat wing appears horizontal on an 
x-y projection. The principal disadvantage is that the first coordinate 
does not appear in two-dimensional eqmtions. The investigation reported 
herein was conducted at Brown University under the sponsorship and with 
the financial assistance of the National Advisory Committee for Aeronautics. 


SUfflOLS 
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velocity potential 
operator 

velocity of sound 
pressure 
density 
entropy 

Cartesian coordinates 
Mach number 
free-stream velocity 
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u,v,w 

a 

y 

Cl 

Cd 

m 

R 

J 

T 

r.BjZ 

t 

(j,01,02 


perturbation velocity con^jonents 
angle of attack or of local slope 
ratio of specific heats - ■ 
lift coefficient 
drag coefficient 

oblique transfonnation parameter • 

hyperbolic distance 

Jacobian 

general homogeneous -wave solution 
cylindrical coordinates 
basic conical coordinate 
conical coordinates iiiithin Mach cone 
conical coordinates outside Mach cone 
arbitrary functions 


“(■'■) jP('t), 7 (t) homogeneous functions 
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Rl 

P 

q 

s 

Co 

ti,t2jS 

T 


fundamental variable 
real part 

rolling angular velocity 
pitching angular velocity 

variable for Jacobi’s, imaginary transformation 
root chord 

geometrical uing parameters 
root thickness ratio 
thickness function . ■ 
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h 

$ 

t 

0 

CO 

k 

h(R) 

n 

r 

J 


velocity potential in unsteady flow 
time 

operator involving solution of wave equation 

frequency 

reduced frequency 

function of hyperbolic distance 

homogeneity parameter 

gamma function 

Bessel function 


ELEMENTS OF LINEARIZED SUPERSONIC FLOW 
Basic Flow Ecjuations 


With the usual assumptions of zero viscosity and heat conduction, 
initial irrotationality and isentropy, uniform stagnation enthalpy and 
fl\rid composition, and no body forces, it may be shown that the fluid is 
always irrotational and isentropic and that there exists a velocity 
potential cp -whose gradient is the vector velocity of the flow field. 
For steady flow this quantity satisfies the equation 


a2^2(p = (vcp)v 



( 1 ) 


■where a is the local speed of sound in the medium, defined as 


a 


2 = 



( 2 ) 


If the flow consists of a uniform flow of velocity V -with velocity 
of sound Uq -which is slightly perturbed -the potential equation may be 
expressed apprcxximately, -wi-bh terms of an order higher "than 0 in the 
perturba-fcion dropped, in the form 


'Pxx + ^jy + (1 “ = 0 


(3) 
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where 


M = V/ao 


(h) 


Two principal assun^jtions are necessary for the validity of this linear- 
ized potential equation. These are: 

(l) The lateral velocity perturbation con^jonents must be small com- 
pared with ao. This requires that the quantity oM be small cou 5 )ared 
with unity idiere a is the local inclination of the flow at any point. 
Thus flows with small values of a are not necessarily linearized for 
large M in the hypersonic range. 


(2) The axial velocity perturbation con 5 )onent must be snail com- 
pared with I V - ao I . This condition would most immediately be expres- 


sible in terms of restrictions ^on the pressure perturbation. However^ 
boundary conditions are usually set in terms of flow inclination, and 
the appropriate condition which appears from the transonic similitude 


theory is that 


a 

|m 2 - Ip/^ 



be small con^ared with unity. 


With such singularities as occur at subsonic leading edges or vortex 
sheet edges one or both of these assumptions fail locally. 

Except in the incompressible or low-speed range M « 1 these two 
assumptions require automatically that the flow inclination a itself 
be small con^ared with rmity, permitting the botmdary conditions to be 
linearized. For a surface locally parallel to the x-axis the boundary 
condition would be 


= V = Va (5) 

where a here is the local inclination of the surface to the principal 
flow direction. With flow in the low-speed range the linearized potential 
equation is valid with nonlinear or iinlinearized boundary conditions. 

The linearized potential eqmtion is considered to be an equation 
for the pertTirbation potential. The three Cartesian derivatives of this 
potential are the three perturbation velocity con 5 )onents u, v, and w. 
The perturbation pressure, or difference between the local pressure and 
that in the free stream, is expressed 


/ u2 . v2\ 
p = _p(vw + —J—J 


(6) 
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The terms in u and t, although nonlinear, are properly included and 
necessary in generalj the skeptic may atteii 5 )t to satisfy the condition 
of zero normal pressure gradient on a circular cone at zero incidence 
1 -iithout them. However, for many cases of ingjortance and, indeed, for 
all those exemplified in this report the additional “planar-system" or 
“mean-surface" assumption holds, and the pressure terms in u and v 
are properly dropped leaving the si np1 i.f ied pressure equation 

p = -pVw (7) 

This additional asstuiption is discussed below. 


Mean-Surface Assumption 

The useful additional restriction here discussed has been previously 
made by most investigators under the terminology of “planar system, “ 
“quasi cylinder," “flat body," “mean surface," and others. In its most 
general form it is assumed that the boundary conditions on a body may 
be satisfied at appropriate points on a mean surface. This mean surface 
must be a part of a general cylinder >rith z-axis directrix; that is, 
must be everywhere parallel to the z-axis. In order for the assunption 
to be valid it is necessary that the body closely approximate the mean 
surface and that the inclination of the body surface to the mean surface 
be small. The case where the body is a wing and the mean surface is the 
xz-plane is the one most commonly encountered. 

The importance of this additional assumption is that it gives a 
direct correspondence between the superposition df solutions T-ri.th the 
same mean surface and the corresponding body shapes. Ihe assumption 
also insures the validity of the simplified pressure equation (7). 


Prandtl-Glauert Transformation 

If with the axial variables unchanged the lateral variables are 
transformed by 


x« 

= - 1 X , 

(8a) 

7’ 

the potential equation becomes 

= \/m 2 - 1 y 

(8b) 

'Px»x« 

'Pyiyl ~ 'Pzz = 0 

(9) 
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which is of the same form as equation (3) with the Mach number equal to 
\[2. Where the planar-system or mean-surface assun 5 )tion is applicable 
it is permissible to set 


cpt 


i/m 2 - 1 cp 


( 10 ) 


this relation leaving the angle of attack, thickness ratio, and so forth 
unchanged under the transformation. However, there must be a change in 
aspect ratio, gap ratio, and so forth 

R« = - 1 R (11) 


and there is a change in the appropriate lift and drag coefficients 

Cl« = \/m2 - 1 Cl (12a) 

Cj)» = \/m 2 - 1 Cj) (12b) 

In the general case where the mean-surface assuit 5 )tion is invalid 
the factor in equation (10) must be (m 2 - l), with the same factor 

appearing in equation (12a), and the factor (M^ - l)^'^^, in equa- 
tion (12 b) . Th e angle of attack and thickness ratio are changed by the 
factor /m^ - 1, 

In the analysis henceforth it is assumed that the Ifech number equals 
/2" so that the primes in equation (9) are dropped. Reduction of the 
results to other values of the Mach nun±)er is carried out by the equa- 
tions given above. 


Oblique Transformation 

The useful oblique transformation is essentially a Lorentz trans- 
formation and was first applied to’ steady supersonic flow by Jones 
(reference 5) . The transformation is in the independent variables 


- (x» - mz*) 

\jj. - m^ 


X 


(13a) 
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2 = ■ ■ - (z» - IJDC*) 

/l - m2 


( 13 b) 


y = y’ 


( 13 c) 


These transformations fona an Abelian group with the inverse trans- 
formation obtained by changing the sign of m. The Cartesian velocity 
components obey the inverse transformation law 

u = - (ui + nM») (Ilia) 

l/l - m2 

1 

w — — — — — (w* + mu') (lijb) 

h - m2 


T = v‘ (1)40) 


Three other properties of this transformation are the invariance of the 
hyperbolic distance, the Jacobian of the transformation being unity, and 
the homographic form of the transformation for x/z: 

R = /z2 - - y2 = - (x»)^ - (y*)^ (l^a) 



(15b) 


(15c) 


\ 
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Symmetry Considerations 

Where- the mean surface is a plane and a planar system exists it is 
convenient to separate the potential with regard to symmetry in the 
y-direction 


<P = 9s 

• (16a) 

<Ps = |[^(y) + 

(16b) 

^a = 

(16c) 


The symmetric part corresponds to the thickness distribution of the. body 
under consideration; the antisymmetric part^ to the mean camber and lift 
distribution. Both qps and cp^ satisfy the potential equation indi- 
vidually and it is possible to treat the two parts separately and 
independently . 

Where the mean surface is symmetrical with respect to the x-axis 
considerations of , lateral symmetry may be applied. With a laterally 
symmetric thickness or camber distribution the potential and pressure 
solutions are laterally symmetric; with a laterally antisymmetric camber 
distribution an antisymmetric lift distribution results. 


Supersonic and Subsonic Edges 

The edges of a mean-surface plan form are conveniently classified 
according to the local-flow orientation. The free-stream velocity may 
be considered to be conposed of two perpendicular velocity con 5 )onents, 
a normal one perpendicular to the edge and a tangential one parallel to 
the edge. The normal one, taken positive from the region off the plan 
form onto ^>he plan form, is divided by the speed of sound in the free 
stream to give a normal Mach riiimber If this quantity is positive 

the edge is termed leading; if negative, trailing. If its magnitude is 
greater than 1, the edge is supersonic; if less than 1, subsonic. If 
Mji = 0 the edge may be called a (subsonic) side edge. 

The nature of the edge as indicated by this ' classification strongly 
affects the nature of solutions for the flow field. The Kutta condition 
preventing an infinite behavior of lifting solutions must be applied on 
all subsonic trailing edges. 
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SOLUTIONS HOMOGENEOUS OF DEGREE 0 


As a preliminary to expressing the generalized conical-flow solution 
the properties of solutions to the wave equation idiich are homogeneous of 
degree 0 are studied. Denoting the quantity of interest by with a 

change to cylindrical coordinates, the wave equation may be expressed 


*^xx ^yy ~ ^zz 

■^rr + J "^r + 4 ° 

r^ 

With the new variable 

t = r/z 

introduced to replace r, the fxmction t must be a function of t 
and 0 alone. The equation for t becomes 


t^(l - t2)T-tt + t(l - 2t2)r^ + T00 = 0 (18) 


In a similar manner the homogeneous solutions to the equation of 
chara ct eristic s 


satisfy the equation of characteristics of equation (18) 

t2(l - t2)T-t^ + T02 = 0 (20) 


The interior and exterior of the Mach cone are treated separately. 

Mach cone interior .— In the case of the Mach cone interior the 
coordinate 
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is introduced and the ■wave and characteristic equations become 

Tgs + T00 = 0 (22a) 

(ts - iT0)(Tg + 1 t 0) = 0 (22b) 

Introducing the con^jlex variable e of the original Busemann linearized 
conical-flow theory 

e = e3q> (s + 16 ) (23) 

■fche solutions to equations (22) may be expressed 

T = f(e) + g(e*) (2iia) 

T = f(e) or g(e*) (2l|b) 

tdiere f and g are arbitrary analytic functions. Of interest for the 
present pxnpose are -the particular functions e, e*, and t*, 

\ihere 




26 
1 + 


( 25 ) 


The principal reason for -the use of this function is that for 0=0 it 
is real and equal to t and for t = 1 it is real and equal to sec 0j 
this gives it the nature of a physical coordinate for a flat wing. 

Mach cone exterior .- In the cage of the Mach cone exierior the 
folio-wing coordinate is Introduced 


CT 


sec 


-It 


( 26 ) 


The wave and characteristic equations -take the form 


■^00 - '^66 = 0 


(27a) 


■^0^ - T02 = 0 


( 27 b) 
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Introducing the characteristic coordinates 

02 ^ = 9 + 0 
02 = 6 - 0 

the solutions may be expressed 

T = f(9i) + g(02) 

. T = f(0i) or g( 02 ) 


(28a) 

(28b) 


( 29 a) 

.(29b) 


The particular solutions of interest are simply 0^^ and 02. 

The Cartesian derivatives of T are functions which are homogeneous 
of degree -1. Functions proportional to these Cartesian derivatives are 
needed which are functions of t alone and hence homogeneous of degree 0 


g(T) _ P(t) _ y(r) 


(30) 


Such functions, provided they exist, are not unique but may be made 
definite by arbitrarily establishing one of them. It is convenient to 
do this by setting 


7(t) = 1 (31) 

In finding these functions it is not necessary to investigate e* or 
because the desired functions will be the complex conjugates of the 
corresponding functions for e and Although a change from g to 

6*, for exaitple, does make a difference in the functions, a change from 
e to T = f(g) vjill make no essential difference, the formulation in 
terms of the different independent variables being, of course, changed. 

Carrying out the indicated Cartesian differentiations yields, after 
some analysis, the following results; 
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T 

a 

P 

7 

6 

1 + 
2e 

. 1 - e2 

^ 2e 

1 

t 

' H I 
1 

r 

Ji- ^ 

1 

01 

-cos 0^ 

-sin 0J 

1 ■ 

02 

-cos 02 r 

-sin 02 



1 


On the Mach cone, e = e^®, ^ = sec 0, 0]^ = 02 = 0, a = -cos 0, and 

p = -sin 0 in all cases. 

GEMEEALTZED CONICAL FLOW 
Fundamental Solutions 


The general solution will be established as a function of the 
"solution poin;^* ^ tei^ of an integral taken over values 

of "the "argument point" (x,y,z) . The designations Tq and r are used 
for the values of the homogeneous quantity of interest t at the solu- 
tion and arguinent points, respectively. A new basic variable is defined 
as a fimction of these poiats Tdiich vanishes when they coincide: 


I 


Xoa(r) + yoP(T) + Zq7(t) 

(32a) 

g = 0 for T = Tq 

(32b) 


■vrtiere the variable vanishes as a result of Euler’s 'formula for homogeneous 
functions of degree 0 


y^y + zTg = 0 (33) 

and equation (3^). The fimdamental idea of Mr. G. N. Ward was the defini- 
' tioh and use of this new basic variable. 
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The generalized conical solution may now be expressed 

«p(3bj7oiZo) g(5,T) dr (3U) 

where the absence of the lower limit indicates that the integral is an 
indefinite one with the usual arbitrary constant. In the exterior of 
the ]^ch cone the general solution is the superposition of solutions of 
the form of equation ( 3 U) with 0p and 82 chosen for t. In the 

interior a real solution is obtained by taking the real part of a solu- 
tion with T equal to e or amounting to a superposition of a 

solution (equation (3U)) with its complex conjugate solution. 

That equation (3U) does give a solution of the wave equation is 
easily ascertained by carrying out the indicated differentiations and 
using the properties of equations (l7a), ( 19 ) j ( 30 ), ( 32 a), and ( 32 b). 
The essential restrictions on g are that it have first and second 
derivatives with respect to I and that it be analytic in x within 
the Each cone. The velocity components are expressed 

^ = 'Pxo = Toxo 


= 9yo = Toy^ g(0,To) 


^O 'Pzo ^OZq S(0,To) 


For most cases of physical interest 

(homogeneous of degree -1) in equations (35) do not appear. 

The solution procedure is divided into three parts: First, the 

solution tpstream from the vertex exterior to the Mach cone, if there 
is any, is to be found, then the solution exterior to the Mach cone 
doimstream from the vertex, and finally the solution within the down- 
stream Mach cone. In each case the solution fomad in the previous step 
provides boundary conditions for the part of the problem at hand. The 
solution in the upstream Mach cone is assaimed null. 


.pair) 


9g(ljT) 


dr (3^a) 


r'^o , , 9g(l,x) 

+ / p(T) 


dr (35b) 


, 9g(l,x) ^ 

+ / y(r) dx 


(35c) 


g(0,x) = 0 and the first terms 
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Oblique Itansformation 

The oblique transformation is applied by equating a known solution 
of the wave equation in the coordinates x, y, and z to a solution in 
the coordinates x*, y*, and z*. This gives a new solution expressed in 
terms of the primed coordinates 


<P( 2 :,y,z) = q)»(x»,y*,z’) 


(36) 


The behavior of the various basic variables under the transformation is 
first obtained. The valuables cos 0, and cos 02 all obey the 

same tiansforraation law as does t cos 0 = x/z in equation (l5c). Thus, 


m 

1 - 


(37a) 


cos 0]_* - m 

cos 0]^ = (37t>) 

1 - m cos 0]^* 


cos 02 * - m 

cos 02 = ' - — (37c) 

1 - m cos 02 * 

The quantity e also follows a homographic transformation, but one 
which is the square root of the other 


e 




m 

1 + /l - m2 


1 - 


m 



e« 


(37d) 


The quantities a and p follow the same transformation laws in all 
cases provided that the relation y = 1 is maintained 

a* + m 


a 


1 + ma’ 


( 38 a) 
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P 


\Jl - 

1 pi 

1 + ma* 


With these results it is easy to show that 


( 38 b) 


I 


\l± - 

! |i 

1 + ma‘ 


( 38 c) 


T^ere the |« is defined in t]ie same way in terms of the primed coordi- 
nates, ■with yi = 1. 

The solution in ■fahe primed coordinates 


<P‘(xoSyo‘ 



gl(|I,Tt) dTl 


(39) 


may now be expressed in tenns of the f ■unction t(t* ) expressed by eq^ua— 
tions (37), and relation (36). The desired relation is given by 


g'(t’,T‘) 


]jl - m2 

1 + ma* 


|i,t(t«) 


dT(T» ) 

dr’ 


(UO) 


Equations (3^) connecting 'the velocity components may be ■used as a check. 


Conical Flow 


The linearized conical-flow theory is ob'tained immediately by 
setting 


g( ^.T) 


df(r) 
d-T . 


(la) 


This leads to the usual expressions for the interior of the Mach cone 
found by Busemann 


(U2a) 


w = Rl |f(e 3 



3F 


MCA TN 2667 


17 


u = El 



1 + €‘ 


2e 



V 




or the analogous expressions in ^ 

w = El [f(e3 



For the exterior of the Mch cone 

w = fi(9i) + f2(02) 


u = - / cos 0x dfi ~ / cos 02 df2 


~ ~J^ sin 0x dfx “ sin 0£ df 2 


Thus the well-known linearized supersonic conical-flow theory 
appears as a special case of the theory here presented. 


Second-Degree Conical Flow 

The case which is of most interest here is that for which the 
potential is homogeneous of degree 2, with 


0l2b) 

0;2c) 

(ll3a) 

(1^3b) 

(Ji3c) 

(i|l4a) 

ihhb) 

ihho) 
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g(ljT) 


|2 df(T) 

‘2 dr 


this corresponding to 


= Rlf(T) 

Noting on the plane y = 0 that 




(U6) 


? = Zq 


^2 = 



(h7) 


the velocity coirponents may be expressed in the interior of the Mach 
cone and on Jq = 0 as 


w = ZqRI! 


}-t) 


df(0 


(U8a) 


u = ZqRI 


df 


V = ZqEI 


-/}(' - 1) 




(U8b) 


(l;8c) 


In the exterior of the Mach cone, again on the plane y^ = 0, the 
analogous expressions are 


w = z. 


u = z. 


= Zr 


7 a - cos 9 i) df^C^i) + J^(l “ "to 8f2(02^ (i; 9 a) 

cos 0i(l - to cos 0^) df;L “ 02^^ ” ^o (hPh) 

- J ' sin 0]^(l - tg cos 0]^) df^ - sin 02(1 ~ to cos 02) df 


(li9c) 


V 
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As is the case with problems in conical flow, the procedure is to 
find a function f so that the prescribed boundary conditions are 
satisfied. The quadratures corresponding to the relations given above 
are needed to con 5 )lete the problem. The exan^jles worked out below serve 
to illustrate the method. 


ROLLING AND PITCHING WINGS 
General Considerations 


Wings of triangular shape in steady rolling and pitching motion are 
now considered, with the apex of the triangle leading and with a straight 
supersonic trailing edge. Only the cases with lateral symmetry are con- 
sidered but cases of yawed wings may be obtained by a suitable use of the 
oblique transfoimation. Such solutions are antisymmetric in the sense 
of the section titled "Symmetry Considerations." For the wing in 
reversed flow the methods of this report do not work. However, the 
pitching- and rolling-moment damping coefficients of the original wing 
apply for the reversed wing, as shown by the reversed-flow theorems of 
Brown (reference 6) and, in still more general form, of Ursell and Ward 
(reference 7) . 

In order for the steady-state equations to be valid the motion must 
be steady in that the local angle of attack at a point on the wing is 
constant. For pitching motion this requires that the flight path be 
curvilinear^ for rolling motion, that the rolling angular velocity be 
constant. In these cases the problem is identical with a problem in 
direct flight with the angle of attack proportional to z or x. 

Where these conditions are not met but where the second derivatives 
with respect to time may be neglected the quasi-steady theory may be 
used. In the cases here disciissed the analysis is similar to that of 
the steady case and the same functions appear. 


Triangular Wings within Mach Cone 

The choice of the function f(D in the cases involving a triangular 
tring with lateral symmetry within the Mach cone is made through the 
following considerations : 

' I 

(1) The function f may have at worst a -3/2-power singularity at 
^ = itj, the leading edges of the. wing ' 

(2) The function df/d^ must have at least a double zero at 5=0 
in order that the integrals for jU and v converge 
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( 3 ) The function f must be bounded for large yalues of 
corresponding to the intersection of the Mach cone with the plane of 
the wing 

(I 4 .) Lateral-sym m etry considerations require that f be an even 
function of ^ in the pitching case and an odd one in the rolling case 


The functions are thus expressed, apart from a multiplicative constant, 
as 


fq(S) 



( 50 ) 


in the pitching case and 

1 ^3 

3 ^ 

2 _ ^ 2 ^/^ 

in the rolling case. The remainder of the problem consists of the 
canning out of the appropriate quadratures in order to evalioate the 
arbitrary constants in terms of the physical problem and to obtain the 
desired velocity component or pressure solutions. 

The quadratures needed to establish the arbitrary constants are the 
appropriate terras from the expression 



(51) 


v(xo,0,Zq) 



or in terms of the pitching and rolling angular velocities 


(52) 



(53a) 


(53b) 
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These Integrals are evaluated by means of Jacobi's imaginary transformation 

i^l - s 2 




iSh) 


in terms of 'which the integrals become 

^1 l/l - s 2 ds 


q = -Ac 


0 


(55a) 


" -W f 


p 

s'^ ds 




where 


(55b) 


(ti«)2 = 1 - ti2 


(56) 


These integrals may be expressed through the classical, although lengthy, 
methods of elliptic Integrals. The resulting expressions are 


Aq = 


-3qti2(l - tq 2 ) 


(1 - 2t3_2)Et(ti) + ti2K'(ti) 


(5?a) 


“3ptp2(i - tj^) 


(2 - ti2)E«(ti) - ti2Ki(tq) 


(57b) 


as 


The pressure on the upper surface p^ is given through equation ( 7 ) 


Pu " 


2 tq2 _ t2 


3ti2(ti2 - t2) 


1/2 


( 58 a) 
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Pu = 


t 

3(ti2 - 


(58b) 


The pitching results may be combined with the well-known conical solu- 
tion to obtain damping coefficients about axes other than that through 
the vertex. The qua si-1 in ear theory may be applied to obtain the appro- 
priate nonsteady results . 

The expressions given here agree with those of reference 8 in all 
pertinent details.. Consideration of the syoometric solutions with n = 2 
with specified pressure on a wing inside the Mch cone is included in 
the appendix. 


Triangular Wings outside Mach Cone 

The triangular wing with supersonic leading edges is characterized . 
by the tangent of the angle which the leading edge makes with the flow 
direction^ suitably reduced according to the Prandtl-Glauert transforma- 
tion and the corresponding value of 6^ 


tq = sec 0Q 

Taking the right-hand side of the wing so that 


( 59 ) 


f 


2 


0 


fq - 9l > Qq 


the solution may be expressed in the form 




2 ^^(2 “ ^ ) 

6(01 - 0o) y 1(00 - e^) 


(tq2 - 1) 


(ti^ - 1) 


3/2 


(60a) 
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fn = p 






(V - 1) 


6(01 - 0 o) - 




(ti2 - 1) 


3/2 


1(00 - 0i) 


(60b) 


where 6 is the Dirac delta function and 1 is the unit step function. 
This leads to the expressions for the pressure on the upper right-hand 
surface of the wing qutside the Mach cone 

Pu = pvq — ~ 3 

(ti2 - 1)3/2 


■tl2 

Pu = P^P ^(2 - tix) 

(ti2 - ir 

for the pitching and rolling cases, respectively. Correspond in g results 
hold for the left-hand side of the wing. 

For the solutions in the interior of the Mach cone it is most con- 
venient to consider the required form for the functions and cpyj- 

in the two cases. A consideration of the singularities leads to the 
results 


( 6 la) 


( 61 b) 


<Pyz 



^ ^1 ^ ^^1 ^ ti \ 

^ ^ - ti ^ C + ^1/ 


( 62 a) 


'Pyx “ ^ 



i - ti 


ti 

c - ti 



(62b) 


whence the form of the function df/d^ may be obtained 
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2h 


VU) - - 3 


t 

r . 



h - 1^ 

^ - tx^ 

(E - 


(63a) 


V(£) - - 1 — 


/l - 


2tn tl tn 

i + +. i. 


it - + t^y 


( 63 b) 


Application of equation (l;8a) then gives the pressure on the i^per sur- 
face of the wing 


Pu = 


- 1 ) 3/2 


+ 


1^(2 - tx^)z - X 

1 -1 1 - txt 

tx + t 

[^(2 - tx^) + £1 

_x ^ 

cos > 

tx + t 1 


(61ia) 


Pu = 


tl2 


2 -n^3/2l- 


ir(ti2 _ 1) 


(z - txx) cos 


tx - t 


(z + txx) 


cos 


-1 


1 + txt 

t^ + t 


(6ijb) 


These results check with those of previous investigators, for exauqple, 
reference 8. It may he readily checked that equations (61|) fit equa- 
tions (6l) at t = 1 and that eqiiations ( 63 ) are consistent with 
equations (60) on the Mach cone. 
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FAMLY OF SWEPT DELTA WINGS 


The methods given above have been applied to the calculation of the 
wave drag of a family of thin wings at zero incidence. These wings are 
symmetric laterally and have a symmetric thickness distributionj thus 
the resulting solution for the velocity potential has the same symmetry 
and there are no lift forces on the wings. The plan-fonn shape is a 
quadrilateral with two of the vertices on the axis, so that the wing 
appears as a delta wing with sweep. 

With the origin at the forward vertex of the wing the leading edge 
is given by the equation 


z 



(65a) 


where t^^ is a parameter indicating the angle between the leading edge 
and the flow direction. If the maximum chord on the central axis is Cq 
and t2 is a similar parameter giving the trailing edge direction, the 
trailing edge is given by the equation 



(65b) 


If the wing has a circiilar-arc profile the upper surface of the 
T-Ting on the plan form must be given by 


Yu = 





t(ix|) 


( 66 ) 


where T is the thickness ratio of the central section and tK^) ~ 1- 
The lower surface of the wing is given sin^jly by the negative of this 
expression. For a wing with constant thickness ratio the function 
must be 

^2 ^1 |x[y^ 

^ 2^1 ' 

Since such a wing cannot be represented by a finite number of homogeneou.' 
flows the choice is made here that 




t =1 


( 68 ) 
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The family of wings 
tion given by 


to be investigated has the thickness distribu- 


= 




(69) 


for which the boundary condition is 


la = 

V 3z 



2z + 


^2 ~ '*^1 
t2ti 



(70) 


It may be noted that such a wing has constant curvature in the axial 
dire c tion j this being 


oVu 

dz^ 


ilT 


(71) 


The solution my be considered to consist of three separate parts 
corresponding to the three terms in the parentheses of equation ( 70 ). 

The first of these is a conical flow of a type much utilized in delta- 
wing calculations and the solution will not be expressed here. The other 
tvro are of the second-degree conical-flow type, one with cppj. = 0 and 


qp- 


yz 


Utv 


and the other with 


qpyz 


= 0 and 


cp. 


yx 


2TV ^2 ~ |x| 

Co t2t]^ K 


( 72 a) 


( 72 b) 


The first of these two solutions of the second degree is, for 
t]^ > 1 and the wing lying outside the Mach cone of the vertex, identi- 
cal with the pitching solution of the preceding section with 


q = 



(73) 
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For the piirpose of coii 5 )ariiig solutions iJith t]^ greater than and less 
than 1 the pressure solution (6Lfa.) may be rewritten as 


l;TpV2 
Pn 


Tr(ti2 - 1)^/^ L 


i-2ti\/ti2 _ 1 


\/z2 - x2 + 


[tiC2 - ti2)z - ^ tan-1 


t]_(2 - tx^)z + X tan“l 


_ 1 /i _ 


1 + tn t 


(7ija) 


For the case where t^ < 1 the corresponding solution is immediately 
obtained from equation (7ija) 


Pu = 


i;TpV2 


Co 


tr(l 


— i — J \/z2 - x2 - 


[7i( 2 - t,2). . .] tann-1 - 


I- 0-1 n \/l - ^ 1 ^ 

[tx(2 - ti2)z + xj tanh ^ ^ 

-J 1 + t^t 


ilkb) 


For the second solution the result has some similarities to the 
rolling case of the preceding section with 


= - ~^2 
Co t]^t2 


( 75 ) 
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the solutions exterior to the Mach cone being identical on the right- 
hand side and of opposite sign on the left. In the interior of the Mach 
cone equation (62b) must be replaced by 





^ ^1 \ 
^2 _ ti2 ^ - ti ^ + t^^y 


( 76 ) 


the new singularity at ^ ~ 0 aj^earing because of the discontinuous 
behavior of the derivative of |x|. Thus equation (63b) is replaced by 


f‘(^) 


£2 

-2ti2 

/i - 

- ti^)' - ti)^ 


h. 

(^ + ti)2 


( 77 ) 


With this formula and relation ( 75 ) the pressure may be expressed for 
t]^ > 1 as 


2 Tpy2 ti - t2 1 


U Co 


V2 n(ti 2 _ 1 ) 3/2 


-2ti|^ti^ - 1 /z2 - x2 + 


ti^(z - tix) tan“t 


_1 \fti2 - l\/l - t2 


1 - tnt 


o 1 A “ 

ti^(z + tpx) tan“t 

1 + tit J 


(78a) 


Correspondingly^ for t]_ < 1 


Pu = 


2TpV2 ti - t2 1 

2tl^l - 

°° *1*2 _ ^^ 2 ) 3/2 

2/ , 1 - ti2 

ti2(z - tix) tanh ^ 

\/l - t2 
tit 

P, . 

ti^(z + t-|x) tanh ^ 

1 + 

\/l - t2 
tit 


- t^ - 


(78b) 
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The remainder of the problem is to multiply these pressures by the 
local angle of attack and to add the results of the integrations of 
these products over the area of the wing. The numerical work involved 
in this procedure which yields the drag is exceedingly long' and tedious 
and no atteirg^t will be made here to repeat any of these calculations. 
The analytical results for the drag of this system of bodies are given 
directly. 

The drag given for this family of bodies by strip theory including 

the \/m^ - 1 factor, which is what the drag should approach for very 
large Mach numbers, may be expressed by 



(79) 


in terms of the drag coefficient. This is less than the two-dimensional 
value for a circular-arc airfoil by a factor of l/2 because the thick- 
ness ratio is not constant, varying linearly from- a zero value at the 
wing tips to the reference value T at the root. For convenience the 
drags are all represented as the ratio of the drag coefficient to this 
reference fictitious value obtained from strip theory. 

The results are given below. The new parameter s is introduced. 


s = ti/t 2 


( 80 ) 


which is geometric, independent of M for a certain wing. For 
1 < ti < 00 and 1 < |t 2 l. 


Cd/CDq 


(1 - 1 

2(1 + 3s2) 1 - 2s - s^ 

^t2 tan“l\ft2^ - 1\ 

n(l - s)^ 1 

1 - s^ s(t2^ - l) 

\ - 1 / 


2s(3 + s2) 1 + 2s - s2 


1 - s2 (ti^ - 1) 


^t]L tan~^[/t]^^ - 
- 1 


-ti 


■ 2s 

tp2 _ 1 


(1 - s 2)2 


s(t2^ - 1) s2(t]^2 _ i)(t2^ _ 


( 81 ) 
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This resTilt must be specialized by suitable limiting processes in several 
special cases: 

Case a: where tp = 1 and t 2 = l/s. 



Case p: vdiere l/t2 = 0 and s = 0, 
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Case (a - 6): where ti = t2 = 1, 

Gj /% = — = 1.11^6 (87) 

o 105n 

These formulas, with the general rever'sed-flow theorem for drag 
(reference l), give the theoretical wave drag con^Jletely for the case 
in which both leading and trailing edges are supersonic. The case in 
■tdiich both edges are subsonic lies outside the scope of the present 
analysis. For the case in which one edge is subsonic -and the other 
supersonic, without loss of generality as far as the drag is concerned, 
the leading edge may be required to be the subsonic one. 

The drag coefficients for the case 0 < t]^ < 1 may be represented 
as the sum of two tenns Gjy^^ + ACp, where Cjq is the drag coefficient 
as calculated by formulas (8l) and (83) with the replacement 



and the ACp is an incremental drag term. This incremental term may 
be expressed for 0 < t]_ < 1 and 1 < |t2| as 
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For the special cases limiti.rig operations are again necessary; 
Case p: where l/t2 = 0 and s = 0, 





2ti(2 - ti^) 


tijpT-ti^ (1 - 1]^) 


3/2 


logf 



By the stratagem of introducing terms of the form | s | the two remaining 
special cases may be expressed together. In these cases the Cp^ term 
is also needed. For 

Case e; where t2 = 1 and t]_ = s, 

and for 

Case k: vrtiere t 2 = -1 and t]^ = -s. 


GDi/CDc 


tr(l - s)^ 


-(!•+ 8s - s^ + 2s3)| 


^ tar±i ^j/l - tj^^ 

i/i - 


■ " — ^1 + Us + s^ + 26s3 - 2s^ 
3 |s| 


(91a) 




° rr/l - s^ |s|(l - s)2L_ 


- i(l - 2s - s2)(l - s)2 


s(^3 - 2s + 8s^ + s3) + 


2s2(1 + s )(2 + 3s 2) logg 


1 - s" 


(91b) 


Although the actual drag coefficient in these last two cases is obtained 
by adding the two terms, it is convenient for con 5 )utational purposes to 
display them separately. 
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Calculations based on the. foregoing forraulas have been made for 

s = ±1, ±\/2/2, ±l/2, ±l/Uj and 0. Figtire 1 shows the geometry of 

the family of wings under consideration. Figm*e 2 shows a plot of 
l/t]_ against l/t 2 to illustrate the way in wMch the parameters are 
related. Figures 3 "to 7 gi"ve the drag curves, plotted logarithmically' 
against t-]_, resulting from the con 5 )utations . The calculated values 

of Cjy'Cp^j including the incremental term for ti < 1 , are given in 
table I . 

It might be noted that in using these formulas for coiT 5 )utation 
the greatest care inust be taken to choose the correct branch of the 
function tan~l. 


APPLICATION TO UNSTEADY-FLOW PROBLEMS 
Basic Eqviations 


The velocity potential in unsteady linearized flow, after the 
application of the Prandtl-Glauert transformation which leads to equa- 
tion ( 9 ), may be expressed 


1 2M / s 

“fxx ^yy - '^zz ^ -“ 3 : '^tt Hz (^ 2 ) 

in which the time coordinate is defined as aQ times the time, so that 
aQ does not appear in the equation. A new fictitious time variable is 
introduced 

T = (m 2 - l)f - Mz ( 93 ) 

which leads to the equation 

+ $yy - < i >22 = i9k) 

The pressinre is expressed by the relations 

. ' ■ = + ‘ ’ (95a) 
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— — “ <5t 
Po^’O 


(95b) 


■vrtiere the derivatives 'Hath respect to z are taken with t constant 
and with t constant, respectively. 


Quasi-Stationary Theory 

The quasi-stationary theory is obtained by assuming that terms 
involving the second derivative with respect to time may be neglected, 
or that only the first two terms in a power series in t need be con- 
sidered. Thus, 

<D = q>o + "bPi (96a) 


® = <Po + ■ 2 «Pl + ■ o '~ (96b) 

m2 - 1 m2 - 1 


Denoting 


q> = o(cpy) (97) 

as the solution to the problem in terms of the specified boundary values 
of cpy' on the plan form of Interest, there results 


$ = 0 





M — — 

“ °(%) " “ C%) 


(98) 


At t = 0 the last term may be dropped but the term in — remains 

m2 - 1 

as a contribution due to unsteady-flow effects. This result also may be 
obtained by taking the first two terms in a frequency expansion of the 
results according to the harmonic theory (cf. the res\ilts of Miles 
(reference 9)). 
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In pitching and rolling motion the solution apart from that due to 
the instantaneous angle of attack may be divided into two parts, one 
due to the appropriate angvilar velocity proper which appears in the 
(po term and one due to the time rate of change of the local angle of 


attack which appears in the -n term. 

- 1 


For steady pitching and 


rolling motion the second term does not appear. This may be checked by 
more accurate analyses which take the rotation of the coordinate system 
into account. 


Harmonic Motion 

The solution is now assTomed to be harmonic of angular frequency co 
in the time. A reduced frequency is introduced 


m2 - 1 

and the new potential function is introduced 

0 « “ e^^*^“^^^cp(x,y,z) 


(99) 


( 100 ) 


This fraction satisfies the equation 


cpjQj. + cpyy - (P 22 - k^q) = 0 (101) 

which reduces to equation (9) for co = 0. The usual boundary condition 
in terms of qPy is of the type 


<Py = (102) 


The pressure is expressed 


P _-icot _ „-iMkz 
Po^o 


(Mcpg - ikcp) 


(103) 



36 


NACA TN 2667 


A soXntion is now sought for equation (lOl) of a form similar to 
that for steady flow 


<p(xo,7o^21o) = KRo)^' dT (lOU) 

1 -Jhere R is the hyperbolic distance 

R = /z2 - - y2 (105) 

and g is considered to be homogeneous of degree n, so that 

^ If '= ng (106) 

os 

A procedure similar to that followed in the steady— state case shows 
that the assumed solution is valid il h satisfies the equation 


2(n +1) o ^ 

hjjj^ + ' hi5 + k^h — 0 


R 




for "vdiich the appropriate solution is 




J i(kR) 




1 

^kRj 


(107) 


(108) 


This approach was suggested by the parallel exposition of Germain and 
Bader (reference 10), 


It is still possible, though the practicality might be questioned, 
to express a solution in terms of a general function of two variables. 
Noting that h may be expressed as 


h = 



(kR sin 0),sin^’*’^ 0 d0 


(109) 


these solutions may be superposed to obtain a single one of the form 
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~ J' J' sin^G^r) dr cos (IcRq sin 0) sin 0 d0 (UO) 


■where the function g is a nexj one suitably constructed from "the old 
ones . 


Brown University 

Protridence, R. I., August I 3 , 1951 
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APPENDIX 

CONSIDEEATION OF ADDITIOML SOLUTIONS OF SYMMETEIC WINGS 


Since the 'writi^ of this report Lomax and Heaslet (reference U) 
have shown that the S3rmmetric solution with n = 2 with specified 
pressure on a i-iing inside the Mach cone is not unique. It is of 
interest to investigate this from the point of view of the theory 
presented herein. The quantity n is equal to ^ + 1 of Lomax and 
Heaslet 's theory. 

In the lift solutions studied in the present investigation, the 
f(^) functions were of the generalized form given by 


dC 



nl 

tn 




where m must satisfy the inequality 


n ^ m ^ 2n - 1 

to prevent imallowable singxilarities on the Mch cone. If m is odd 
the solution is laterally symmetric j if m is even the solution is 
laterally antisymmetric. The number of solutions available is in each 
case equal to the number needed to cover the possible behavior of cpy 

on the wing, so that no eigensolutions are indicated. 

In the symmetric or thickness case the solutions are of the form 


ifM 
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Inhere m now satisfies the ineqxiality 

n ^ IB ^ 2n 

In this case for any n there is one solution available in excess of 
the nun±)er needed to cover the possible behavior of cp^ on the -wing, 
thus yielding one eigensolution for each value of n. For n odd the 
extra solution appears with m odd and is correspondingly laterally 
antisymmetric. For n even the extra solution appears with m even 
and is laterally symmetric. In the case n = 2 the nonuniqueness is 
that foiond by Lomax and Heaslet and the theoiy of this report provides 
a complete check of their result. The characteristic solution still 
exists without the lateral symmetry as long as both edges are subsonic 
leading edges. 

In the conical case n = 1 the characteristic solution is laterally 
antisymmetric and has not, to the authors’ knowledge, been previously 
found. For this case the solution is expressed by 



giving the normal velocity distribution 



This gives for the half -thickness distribution of this wing 

u-l ^■*'0 

y = X cosh 

X 

As was to be expected this distribution of thickness is negative on one 
side. It may be superposed upon the known solution of constant pressure 



to give a family of conical wings satisfying the same boundary conditions. 
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‘ CALCULATED VALUES OF 



a» 

>*1 


TABLE I 

Cj/Cp^ OR Gjy^jCiy^ + AC^yCi)^ 


\J2 1 


^/2 1/2 \/2/k l^k 


1.1156 0.8809 0.7i;3l O.5097 

1.1101 1.1228 .8711; O.6I3I 

1.1061; 1 . 155 k .7287 

1.1073- .861;3 


1.1061; 1 . 155 k 

1.1073- 

1.1630 

1.1138 


1.1285 1.2122 .86i;5 .7188 .5989 .l;9l;0 


1.1552 . 

1.3008 

1.2008 .8771 

1.2a? l.kl 5 k 

1.10l;2 1.8035 .9ll;7 


.7168 


.9921; I .7301 


.5898 
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FIGURE I.-WING GEOMETRY. 
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